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Abstract 

We use a microscopic three-cluster model with the a + a + n configuration to analyze with 
astrophysical interests the photodisintegration cross section of 9 Be. The valence neutron in 9 Be is 
treated in the molecular orbit model (MO), including continuum states by the box discretization 
method. It is shown that good agreement is achieved with the recently measured B(E1) transition 
probability from the ground state to the first l/2 + state when the box size is enlarged, resulting 
in reasonable reproduction of the experimental data for the photodisintegration cross section. The 
Coulomb dissociation of 9 Be is also discussed. 
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The 9 Be nucleus has been popular especially for the past few years. Firstly, it has been 
often used as a target for projectile fragmentation reactions which involve the radioactive 
beams [|IJ. In order to analyze the experimental data, a new neutron- 9 Be optical potential 
was also constructed M. Secondly, the nucleus has been utilized to investigate the role of 
break-up process in heavy-ion subbarrier fusion reactions HQ- It was reported in Ref. 
H that the complete fusion cross sections for the 9 Be + 208 Pb reaction were hindered by 
68% at energies above the Coulomb barrier compared with those expected from a simple 
one dimensional potential model. This large suppression of complete fusion cross sections 
was attributed to the low energy thresholds for breakup of 9 Be into charged fragments, i.e., 
Q = -1.57 MeV for 9 Be -> a + a+n and Q = -2.47 MeV for 9 Be -> a+ 5 He. The subbarrier 
breakup reaction of 9 Be was also measured subsequently, observing large a particle cross 
sections even at energies well below the Coulomb barrier 

9 Be is interesting also from the nuclear structure point of view |5|, |, [7|, ||, ||. The 
prominent a cluster structure is well established for this nucleus, and furthermore, it is a 
loosely bound nucleus in a sense that none of the two-body subsystems (a-a or a-n) is 
particle stable. These properties make 9 Be a good testing ground for any nuclear structure 
model applying to heavier Be isotopes. Based on a similar idea as the linear combination of 
atomic orbitals (LCAO) for hydrogen molecules in quantum chemistry, the molecular orbit 
model (MO) was developed for the 9 Be nucleusf^J. This model has recently been applied 
with success to other heavier isotopes as well |J [ID], [111 [TIJ. For instance, it was pointed out 
that the presence of valence neutrons enhances the a cluster structure of 10 Be by occupying 
the 1/2+ (er-) orbit in the sd-shell Q|. 

Another important aspect of the 9 Be nucleus is its role in the r-process nucleosynthesis for 
heavier elements in the neutrino-driven wind formed in core-collapse supernovae. In such 
an exploding environment with abundant neutrons, the three-body reaction a(an, 7) 9 Be, 
followed by 9 Be(a,n) 12 C, plays a key role to bridge the unstable mass gaps at A=5 and 



8 compared with the triple-a reaction |T3|, [TJJ]. The reaction rate for the former process 
is governed by the low energy resonances of 9 Be, which can be determined experimentally 
by measuring the cross section for the inverse reaction, i.e., the photodisintegration process 
9 Be(7,an)a WM. Utsunomiya et al. recently performed such measurement using laser- 
induced Compton backscattered 7 rays JT5|] . They obtained the B(E1) transition probability 
from the ground 3/2~ state to the first l/2 + state (the first excited state) almost twice as 
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large as the previous measurement which used inelastic electron scattering, and the resonance 
energy of the first l/2 + state which is slightly shifted upperward by 0.064 MeV |L5 . 



From the theoretical point of view, the treatment of the l/2 + state in 9 Be has been a 
challenge, since it is not bound but lies above the 8 Be+n threshold only by 19 keV (E x = 1.68 
MeV). The aim of this paper is to present a detailed analysis of this state based on the 
molecular orbit method, and apply it to the photodisintegration of 9 Be in order to analyze 
the new experimental data of Utsunomiya et al. Although Descouvemont has done similar 
studies based on a microscopic three-cluster model and obtained nice agreement with the 



experimental data for the B(E1) value as well as the photodisintegration cross sections |16 
we particularly focus on the dependence of the B(E1) value on the properties of the wave 
function for the valence neutron. Therefore, to some extent, our study is complementary to 



Ref. |fL6[ . In addition, we also give a brief discussion on the subbarrier Coulomb dissociation 
of 9 Be. 

We introduce a microscopic a+a+n model for 9 Be. We first consider an intrinsic wave 
function as a basis where two alpha particles and the valence neutron are located at ±d/2 
and R n against the origin, respectively. The total intrinsic wave function is fully antisym- 
metrized, and the single particle wave function for all the nucleons is expressed by a Gaussian 
function. Each a cluster contains four nucleons, and its wave function reads 

<P {a) {r ph f pl ,f nh r nl ) 

= G rS^) G rS^pI) liS^)G R a {rnl)Xp^XplXn^Xnl, (1) 

where x denotes the spin-isospin wave function. Here, G represents a Gaussian function, 

G A& = (^7i ex P[-( f - i? ) 2 / 2 ^], (2) 

and R a = ±d/2. The wave function for the valence neutron is also expressed by a Gaussian 
(§) with R = R n . The total intrinsic wave function is then projected onto the eigen-states 
of angular momentum J and parity, which we perform numerically. Note that for 9 Be the 
MO picture can be introduced by merely performing the parity projection, since the a-a 
core is always an even-parity system. 

The spectra of 9 Be are obtained by superposing the basis states with different a-a dis- 
tances d as well as neutron distances R n . The coefficients for the superposition are deter- 
mined by the variational principle. To this end, we descretize R n , and include from zero 



to |-R™ ax | = 15fm with a step size of typically 2 fm. The a-a distances d and the relative 
angle between d and R n are also discretized similarly. For the oscillator parameter in the 
Gaussian function we use s = 1.46 fm in the following calculations. The adopted effec- 
tive nucleon-nucleon interaction is the Volkov No. 2 with the exchange parameters M = 0.6 
(W = 0.4) for the central part, and the G3RS spin-orbit term for the spin-orbit part as in 
Ref. 0. All the parameters of this interaction were determined from the a + n and a + a 
scattering phase shifts. 

Our interest in this paper is to compute the photo reaction cross sections of 9 Be. For this 
purpose, we assume that the 9 Be eventually breaks uptoa + a + n with 100% of probability 
once it is excited above the threshold for neutron emission. With this assumption, the 
photodisintegration cross sections are equivalent to the photo absorption cross sections. 
Since the calculations are based on the bound state approximation, the continuum states 
appear as discrete states. For unpolarized photons, using the Fermi's Golden rule, the El 
photoabsorption probability per unit time from a bound state i to / reads W? 



wSJ(E 7 ) = ^ ~ ^\(f\m 2 6(E f - Ei - E,), (3) 
where E 1 is the photon energy, and Ei and Ef is the energy of the state % and /, respectively. 

/■pi \ 

The photo absorption cross sections are obtained by dividing the probability w^f by the 
photon flux, c/(27r) 3 , 

aSj(E,) = ^f^ B(E1; i - /) 5(E f - E t - E,). (4) 

Here we have carried out the summation for different m substates for the states % and /. 
In order to mimic the neutron emission from the excited states of 9 Be, we smear the delta 
function in Eq. (f|) by the Breit-Wigner function as 

(Eij/px 167r 3 E*y . . r/2vr 

The total cross sections are then obtained by summing up all the final states /. 

In the present model, the energy of the ground 3/2~ state of 9 Be and the ground + state 
of 8 Be is calculated as —56.3 MeV and —54.7 MeV, respectively. Since it is experimentally 
known that the 9 Be nucleus is bound from the 8 Be+n threshold by 1.7 MeV, the present 
model provides the correct binding energy for the ground state of 9 Be within 0.1 MeV 
accuracy, which is remarkably good. The convergence feature of the energy of the ground 



TABLE I: The energy convergence of the ground 3/2~ state and the first l/2 + state of 9 Be. The 
Gaussian center of the valence neutron is incorporated up to |i?™ ax | = (I): 2.5 fm, (II): 6.0 fm, (III): 
10.0 fm, and (IV): 15.0 fm, measured from the center of mass of the a-a core, respectively. The 
energy of the + ground state of the a-a core is calculated to be —54.7 MeV. The B(E1) value 
from the ground 3/2~ state to the first l/2 + state is also listed for each model space. 





\R™ X \ (fm) 


3/2' (MeV) 


1/2+ (MeV) 


B(E1) (e 2 fm 2 )| 


(I) 


2.5 


-54.5 


-50.3 


0.022 


(II) 


6.0 


-55.8 


-53.1 


0.032 


(III) 


10.0 


-56.3 


-54.0 


0.042 


(IV) 


15.0 


-56.3 


-54.3 


0.048 



state as well as the unbound first l/2 + state (just above the 8 Be+7i threshold) with respect to 
the maximum spatial separation |-R™ ax | for the valence neutron is listed in Table I, together 
with the corresponding B(E1) values. Here, the Gaussian center of the valence neutron is 
incorporated up to (I): 2.5 fm, (II): 6.0 fm, (III): 10.0 fm, and (IV): 15.0 fm, measured from 
the center of mass of the a-a core, respectively. Since the ground 3/2~ state is bound from 
the threshold by only 1.7 MeV and it thus has a rather long tail, we see that at least |-R™ ax | 
= 10.0 fm is needed for the quantitative description of the ground state of 9 Be. Note that 
|i?™ ax | = 4 fm was employed in the previous MO calculation for 9 Be 0. For the unbound 
l/2 + state, the main component is the s-wave, which is free from the centrifugal barrier, 
and an extended model space is required. We find that the model space (IV), where the 
wave function of the valence neutron is solved up to 15 fm from the core, is needed to reach 
the convergence. 

For the a + a + n system, the El transition is allowed only as a consequence of the recoil 
effect of the a-a core caused by the valence neutron around it, since we do not introduce the 
effective charge. Notice that the B(E1) value is exactly zero for the a-a core without the 
valence neutron. Therefore, it is crucial to take a large enough model space for the valence 
neutron in order to obtain a quantitative prediction for the B(E1) value. The calculated 
B(E1) value increases by more than a factor of two when the model space is enlarged from 
(I) to (IV), leading to a value which is consistent with the recent experiment by Utsunomiya 
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TABLE II: The excitation energy and the El transition probability from the ground state to the 
excited states used to obtain the photodisintegration cross sections. The (*) denotes the resonance 
states. Only those states whose excitation energy is below 5 MeV are shown. 



r 


E x (MeV) 


B(E1)| (e 2 fm 2 ) 


1/2+ (*) 


2.05 


0.048 


1/2+ 


3.55 


0.0062 


3/2+ (*) 


4.42 


0.0088 


3/2+ 


4.91 


0.018 


5/2+ (*) 


3.13 


0.028 


5/2+ 


4.93 


0.035 



et al. (0.054 ± 0.004 e 2 fm 2 ) ||. 

Adopting the B(E1) value for the model space (IV), we show in Fig. [I] the photodisin- 
tegration cross sections of 9 Be as a function of the photon energy E~. The experimental 



data are taken from Ref. [15|]. Table II lists the states included in the calculation, which 
are obtained by the present MO method. For the width V in eq.([5]), we used the experi- 
mental values for the first 3/2+ and 5/2+ resonance states, while we arbitrarily chose V = 
1 MeV for non-resonant continuum states. For the width for the first 1/2+ state, we use 
the energy dependent width given in Ref. \TE\. The dashed line shows the results of the 
present calculation thus obtained. It appears that the low energy cross sections are sensitive 
to the energy of the lowest 1/2+ state as well as the B(E1) value to this state. We found 
that the energy of the 1/2+ state was rather sensitive to the box size, even after the B(E1) 
value was converged. Since we could not obtain the converged energy for the 1/2+ state 
with controlled numerical accuracy, we instead used the experimental value for it (E x =1.78 
MeV), although we used the calculated value for the B(E1) (see the solid line). Besides 
this point, we see that the agreement of the theoretical results with the experimental cross 
sections is satisfactory. 

We have also applied the present calculation to the Coulomb excitations of 9 Be. In Ref. 
||, it was shown that, although the energy dependence was well reproduced, the measured 
probability for the Coulomb excitation for the 9 Be + 208 Pb reaction were underestimated 
by a factor of two by a calculation which included the first 1/2+ and 5/2~ states together 
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FIG. 1: The El photo-disintegration cross sections of 9 Be as a function of photon energy obtained 



in the molecular orbital model (the dashed line). The experimental data are taken from Ref. [15]. 
The solid line is the result of the molecular orbital model but with use of the experimental energy 
for the lowest l/2 + state. 



with the experimental values for the excitation energies and the electromagnetic transition 
probabilities. The underestimation of the Coulomb excitation probability is likely attributed 
to the lack of non-resonant continuum states in the calculation. We therefore repeated 
the calculation by including the states shown in Table II. It was found, however, that the 
number of continuum states which we obtained by the present MO calculation was too small 
to resolve the discrepancy. This is certainly a limitation of the box discretization, and one 
would need the proper scattering boundary condition for the continuum states in order to 
reconcile the discrepancy. 

In summary, the photodisintegration cross section of 9 Be newly measured by Utsunomiya 
et al. was successfully reproduced by a microscopic a-cluster model. It was found that the 
B(E1) transition from the ground state to the first l/2 + state, which plays an important role 
in explosive nucleosynthesis, has a strong dependence on the model space for the valence 
neutron. Incorporating the neutron orbit up to around 15 fm from the 8 Be core, we obtained 
the B(E1) value consistent with the new experimental value, which is almost twice the 
previous measurement. 

We thank H. Utsunomiya for useful discussions and providing us with the experimental 
data in a numerical form. K.H. thanks D.J. Hinde and M. Dasgupta for useful discussions. 
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